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Abstract 

We conjecture that the free-fermion part of the eigenspectrum ob- 
served recently for the SUq{N) Perk-Schultz spin chain Hamiltonian in 
a finite lattice with q = exp(iTr{N — 1)/N) is a consequence of the exis- 
tence of a special simple eigenvalue for the transfer matrix of the aux- 
iliary inhomogeneous SUq{N — 1) vertex model which appears in the 
nested Bethe ansatz approach. We prove that this conjecture is valid 
for the case of the SU{3) spin chain with periodic boundary condition. 
In this case we obtain a formula for the components of the eigenvec- 
tor of the auxiliary inhomogeneous 6-vertex model {q = exp(2i7r/3)), 
which permit us to find one by one all components of this eigenvector 
and consequently to find the eigenvectors of the free-fermion part of 
the eigenspectrum of the SU (3) spin chain. Similarly as in the known 
case of the SUq{2) case at q = exp(i27r/3) our numerical and ana- 
lytical studies induce some conjectures for special rates of correlation 
functions. 



1. Introduction 

It was found very recently that part of the eigenspectrum of some quan- 
tum spin models is given by a sum of free-fermion quasienergies [|l H]. In 
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particular, the ground state energy of the SUq{N) invariant Perk-Schultz 
Hamiltonian 



L-l 



where 



N-l N-l 



a=0 h=a+l 

and E"-^ are N x N matrices with elements {E"'^)cd = S^S'^, is given by 

, TT , sin 7r(2ri + 1)/2L , , 

B„= 1 + 2(1-1 + „)oos(-) i^^^^ (2) 

for the special value of the deformation parameter 

/in(N-l)\ 
. = exp(^^). (3) 

The lattice size defining the quantum chain is L and the parameter n in 
is given by the integer part of L/N. 

Moreover an amazingly simple formula was found (conjectured) for the 
ground state energy of SUq{N — 1) Hamiltonian for the same value of defor- 
mation parameter 

Eo = -2(L-l)cos(^). (4) 

While for the SUq{N) model many eigenenergies can be described as 
the sum of free-fermion quasienergies, in the case of the SUq{N — 1) there 
exists a single energy level in this class, namely the ground state energy 
given in (^). Actually this state can be considered as a special state of 
the SUq{N) model in the sector where only (A^ — 1) classes of particles are 
present (see conjecture 3 of [@]). We intend to show in this paper that the 
corresponding wave function possess nice combinatorial properties and its 
components play a very important role in the nested Bethe ansatz approach 
being used for the construction of all other free-fermion eigenstates. This 
comes from the fact that using the nested bethe ansatz (NBA) method for 
the SUq{N) invariant Hamiltonian (|lD we obtain (see, for example, [Q) an 
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auxiliary transfer matrix of the inhomogeneous SUq{N — 1) invariant vertex 
model. We conjecture that this matrix has an unique factorizable eigenvalue 
that reduces some of the NBA equations into a very simple form leading to 
the free-fermion-like structure. As a result we obtain the free-fermion part 
of the eigenspectrum of the Hamiltonian (|l]). All related eigenvectors can be 
found from the knowledge of this unique eigenvector of the transfer matrix of 
the inhomogeneous SUq{N — 1) vertex model. In the homogeneous case the 
transfer matrix commutes with the SUq{N — 1) invariant Hamiltonian and 
this special eigenvector reduces to the ground sate eigenvector with energy 

As far as a periodic boundary case is concerned, we have found the 
free-fermion spectrum only for the SU{3) Perk-Schultz model with q = 
exp(2z7r/3) (we do not consider here the free-fermion point of the SU{2) 
model, or equivalently the standard XY spin chain). We show below that 
the free-fermion part of the eigenspectrum of the SU{3) model with periodic 
boundary condition can be explained by the existence of an unique factoriz- 
able eigenvalue for the transfer matrix of the inhomogeneous 6-vertex model 
with q = exp(2i7r/3). We prove the existence of this eigenvalue using different 
methods described in papers 0,0. Comparison of these approaches sheads 
a new light on recently discovered combinatorial properties of the ground 
state wave function for the odd length XXZ spin model with (A = —1/2) 
0) [0)101 (see also [§] for similar results for the XYZ spin chain). 

2. The SUg{N) invariant model 

The Hamiltonian (|l|) describes the dynamics of a system containing 
classes of particles (0,1,. . . ,N-1) with on-site hard-core exclusion. The num- 
ber of particles of each specie is conserved. Consequently we can separate 
the Hilbert space into block disjoint sectors labeled by {no,ni, . . . ,n^_i), 
where rii = 0,1,..., L is the number of particle of specie i (i=0,l,...,N-l). 
The Hamiltonian has a Sn symmetry due to its invariance under the per- 
mutation of particles species, that imply that all the energies can be ob- 
tained from the sectors {no,ni, ... ,11^-1), where no < ni < ■■■ < un-i 
and uq + 111 + . . . + n^-i = L. Moreover the quantum SU{N)g sym- 
metry of H implies that all energies in the sector {nQ,n[, . . . ,n']^_i) with 
^0 — ''^'1 — ■ ■ ■ — ""-Af-i ^'^e degenerated with the energies belonging to 
the sectors {no,ni, ...,nN^i) with rio < < • • ■ < n^^i, if < uq and 
+ n[ < no + rii and so on up to + ri'^ + ■ ■ ■ + n'^_2 < no + ni + - ■ ■ + nN-2- 
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The nested Bethe ansatz equations (NBAE) for the SUq{N) Perk-Schultz 
model (|1|) are given by (see e. g. ||, [1^, [n|, |[) 

Pk Pk-l Pk+1 

j=l,jj^i j=l j=l 

where k = 0, 1, . . . , iV — 2 and i = 1,2,... The integer parameters pk 
depend on the particle occupation numbers {rii}: 

k 

Pk = ^ni, /c = 0, 1, . . . , - 2, p-i = 0, pn-i = L, (6) 

i=0 

and the functions F{x,y) and f{x,y) are defined by 

cos(2y) - cos(2x - 2r]) cos(2y) - cos(2x - rj) 

F{x,y) = — — — — — , f{x,y) = — — — ■ — -. 7 

cos(2y) — cos(2x + zr]) cos(2y) — cos(2x + rj) 

In the NBAE (^ we have variables of different classes. The number of vari- 
ables uf'^ of class k is equal to pk and the variables u[^~^^ = (i = 1, . . . , L). 
The whole system of NBAE consists of subsets of equations labelled by k 
and containing precisely pk equations (fc = 0, 1, . . . , A^ — 2). 

The eigenenergies of the Hamiltonian (|Tp in the sector (^o, rii, . . . , ni\f_i) 
are given by 

^_ y/^ q ^1 sin(Mj-V2) sm{uj+ri/2) \ 
t:^V a sm(ui + n/2) sm(ui-n/2)J 



where to simplify the notation p = Pn-2 and Uj = u^^~'^\j = 1,2, ... ,p. 

All the solutions of NBAE (^) described in this and our previous paper 
satisfy the additional "free-fermion" conditions (FFC) 

/^(m„0) = 1, 1 = 1,..., p. (9) 

Consequently from (§) and (H) the corresponding eigenenergies of the Hamil- 
tonian (|I]) are given by 



E = -2 J2 i^-cosT] + cos^j, 1 < kj < L -1 (10) 



P^'V nk, 

I — cos T] + COS 

where {kj} is any set of distinct integers in the range 1 < kj < L — 1. 
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In the derivation of the NBAE (|^) it is necessary to find the eigenvalues 
of an auxihary matrix that corresponds to the transfer matrix of an inhomo- 
geneous SUq{N — 1) invariant vertex model on the square lattice of width 
p in the horizontal direction. These eigenvalues enter into the NBAE (see 
equation (51) from 0)0 

= r\u„0) n F{u„u,), z = l,...,p. (11) 

We see that the FFC @ could be explained if there exists an special eigen- 
value for this transfer matrix: 

aSj^'\u) = F-\u,u) flFiu,u,). (12) 

This observation can be formulated in the following conjecture. 

Conjecture: let q = exp and consider the inhomogeneous SUq{N— 

1) invariant vertex model on the square lattice with p columns, where p = 
{N — 1) k + r, and < r < N — 2. The inhomogeneity of the model in the 
horizontal direction are fixed by the vertical rapidities Uj, j = 1, . . . ,p. The 
row-to-row transfer matrix, depending on the spectral parameter u (horizon- 
tal rapidity), has a special factorizable eigenvalue given by formula (|7^j. The 
corresponding eigenvector belongs to the sector S = {no, ni, . . . , nN-2} where 
Ui = k, i = 0, . . . ,r — 1 and = k -\- 1 for i = r, . . . , N — 2. 

3. The SU{3) Perk-Schultz model with periodic boundary 

The SU (3) Perk-Schultz model |]12| is the anisotropic version of the SU (3) 
Sutherland model [0], with the Hamiltonian, in a periodic L-site chain, given 
by 

L 

Hg = J2 (^L,L+1 = Hl,i), (13) 

1 2 

Hij = - E E (^f 4" + E^'^Ef - qE'^'^Ef - l/qEfEf). 

a=0 h=a+l 



^AU details related to the construction of the transfer matrix can be found in reference 
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In our previous paper (see |]T|] for details) it was shown for the periodic 
model that: 

The Hamiltonian ( |7^ with L sites at q = exp(2z7r/3) has eigenvectors (not 
all of them) with energy and momentum given by 

E, = -E(l + 2cos^), (14) 

Pi = tT.J (15) 

with I being a subset of X unequal elements of the set {1,2, The 
number I has to be odd I = 2k + 1 and the sector of appearance of the above 
levels is Sk = ik,k + l,L-2k-l),0<k<{L- l)/2. 

The corresponding solutions of the NBAE were described in p| and we 
intend now to present a procedure that allow us to calculate the related wave 
functions. This procedure is based on the coordinate Bethe ansatz method 
and we follow here reference ||14|| . 

Due to the conservation of particles the total numbers of particles uq, rii 
and n2 = L — uq — rii on classes 0, 1 and 2 are good quantum numbers, 
and consequently we can split the associated Hilbert space into block dis- 
joint sectors labeled by the numbers no and ni. We consider the eigenvalue 
equation 



H\nQ,ni >= E\no,ni > (16) 

where 

|no,ni >= ^'^f{xi,Qi; . . . ■,Xn,Qn)\xi,Qi; ...■,Xn,Qn> (17) 

{Q} {^} 

and n = no + ni. Here \xi, Qi, . . . ; x„, Qn > means the configuration where 
a particle of class Qi {Qi = 0, 1) is at position L). The 

summation {Q} = {Qi, . . . , Qn} extends over all (0,1) sequences in which uq 
terms are and ui terms are 1. The summation {x} = {xi, . . . ,Xn} runs, e 
n increasing positive integers with 1 < xi < ■ ■ ■ < x„ < L. Before getting 
the results for general values of n let us consider initially the simple cases 
where we have 1 or 2 particles. 
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n=l. For one particle on the chain (class or 1), as a consequence of 
the translational invariance of (|13D it is simple to verify directly that the 
eigenfunctions are the momentum-/c eigenfunctions 

L L 

|1,0>=E/(^'0)I^'0> |0,l>=^/(x,l)|a;,l> (18) 



x=l x=l 

with 

/(x,0) = /(x,l) = e^'=^ ^=T^' ^ = 0,l,...,i^-l (19) 
and energy given by 

E = e{k) = -e'^ -e-'^ + q + l/q. (20) 

n =2. For two particles of classes Qi and Q2 {Qi,Q2 = 0,1) on the 
lattice, the eigenvalue equation (|^) gives us two distinct relations depending 
on the relative location of the particles. The first relation comes from the 
amplitudes where a particle of class Qi is at position xi and a particle Q2 is 
at position X2, where X2 > Xi + 1. We obtain in this case the relation 

Ef{xi,Qi]X2,Q2) = - f{xi-l,Qi;x2,Q2) - fixi,Qi]X2 + l,Q2) 

- f{xi + 1, Qi] X2, Q2) - f{xi, Qi; X2 - 1, Q2) 
+ 2(g + l/g)/(xi,Qi;x2,Q2). (21) 

This last equation can be solved promptly by the ansatz 

/(xi,gi;a:2,g2) = e*'^^^e*'=^^^ (22) 

with energy 

E = e{k,) + e{k2). (23) 

Since this relation is symmetric under the interchange of ki and k2, we can 
write a more general solution of (21) as 

fix^,Q^■,X2,Q2) = ^ A^-g^e^('=^i^i+'=^2^^) 

p 

= ^^l''32gi(fcia;i+fc2a;2) _j_ j^QuQ2 ^i{k2Xi+kiX2) (24) 
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with the same energy as in (^). In the last equation the summation is over 
the permutations P = Pi,P2 of (1,2). The second relation comes from the 
amplitude where X2 = xi + 1 (matching condition). In this case instead of 
([2T| ) we have 

Efixi, Qi, xi + 1, Q2) = -f{xi - 1, Qi\ xi + 1, Q2) - f{xi, Qi, xi + 2, Q2) 
-fixi, Q2; Xi + 1, Qi) + {2q + l/q)fixi, Qi, Xi + 1, Q2) Qi < Q2 

Ef{xi, Qi, xi + 1, Q2) = -f{xi - 1, Qi; xi + 1, Q2) - f{xi, Qi, xi + 2, Q2) 
+ iq + l/q)fixi,Qr,Xi + l,Q2) Qi=Q2 (25) 

Ef{xi, Qi, xi + 1, Q2) = -f{xi - 1, Qi\ xi + 1, Q2) - f{xi, Qi, xi + 2, Q2) 
-f{xi, Q2; xi + 1, Qi) + {q + 2/q)f{xi, Qi, Xi + 1, Q2) Qi > Q2- 



Since the ansatz (^^ with (^) is also a solution of (pT]) with X2 = xi + 1, 
we obtain from (^) 

/(xi, Qi, xi, Q2) + f{xi + 1, Qi\ xi + 1, Q2) = 

^/Qf{xi, Qi, xi + 1, Q2) + f{xi, Q2] xi + 1, Qi) Qi < Q2 

f{xi, Qi, xi, Q2) + f{xi + 1, Qi\ xi + 1, Q2) = 
(g+l/g)/(xi,Qi;xi + l,Q2) Qi = Q2 (26) 

/(xi, Qi; xi, Q2) + f{xi + 1, Qi\ xi + 1, Q2) = 

q f{xi, Qi; xi + 1, Q2) + /(xi, Q2; xi + 1, Qi) Qi > Q2. 

If we now substitute the ansatz ( P^ into these equations the constants A^j'*^^ 
and y4^^'*^^, initially arbitrary, should now satisfy 

p 

E^a,P2^p5, = g = l,2 (27) 
p 

p 

where 

ap^ P^ = 1 + e^'^^i+^^^a - (g + g-^) e^'^^^ . (28) 



8 



The system (^) consists of 3 equations. The second equation can be easily 
rewritten as 



aQ,Q 



^P2,Pl 



^Pl,P2 

and combining the first and the third equations we obtain 



A 



0,1 

Pl,P2 



-A 



P2,Pl 



^Pl,P2 



+ 
+ 



D«fcPl _ g«'=P2 



(^Pl,P2 



-A 



1,0 



-A 



^Pl,P2 



0,1 

P2,Pl 



;i_^eifcpi)(l_g-ie»fcp2) 

T ^P2 ,Pi ■ 

(^Pl,P2 



The equations (p9|-pl|) can be written in a compact form 

Afj: = - E s^i:2iikp,,kpM%'S (Qi,Q2 = o,i) 

q;,q^=o 

where we have introduced the S-matrix 



^0:0(^1,^2) 

So,i{ki,k2) 
Soj{ki,k2) 



0"2,1 



S{i{ki,k2) 

0"l,2 

(1 - ge^'^^)^! _ g-igifci) 

Cri,2 

(1 - ge^'^Oll - q'^e'^^] 



0"1,2 



'S'l'o(^l5 ^2 



(^^ik2 _ ^iki^ 
0"1,2 



(29) 



(30) 



(31) 



(32) 



(33) 



The equations (|29| - |3T|) or ( p2D do not fix the "wave numbers" ki,k2- In 
general, these numbers are fixed due to the cyclic boundary conditions: 



/(xi, Qi] X2, Q2) = f{x2, Q2] xi + L, Qi 



(34) 
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which from (^) give us the relations 

aQi,Q2 aQ2,Qi Akp. L (QK\ 

^Pi,P2 ~ ^P2,Pi ^ ■ y"^"^) 

This last equation, when solved by exploiting ([3^ ) and (0), gives us the 
possible values of ki and /c2, and from ( P5D the eigenenergies in the sectors 
containing 2 particles. Instead of solving these equations for the particular 
case 77. = 2 let us now consider the case of general n. 

General n. The above calculation can be generalized for arbitrary values 
of no and rii of particles of classes and 1, respectively {tiq + rii = n). The 
ansatz for the wave function ( ]T7| ) becomes 

/(xi, Qi; . . . ; Xn, Qn) = E A^^\:;;^;e^(^^i^i+ - +'=^"^") (36) 

p 

where the sum extends over all permutations P of the integers 1, 2, . . . , n. It is 
simple to see that the relations coming from the eigenvalue equation (^) for 
the components Qi; • • • ; Qn > where Xj+i — Xj > 1 for i = 1, 2, . . . , n 
are satisfied by the ansatz (^Gj) with energy 

n 

E = J2e{k,). (37) 

On the other hand if a pair of particles belonging to classes Qi, Qi+i is located 
at positions Xi, Xj+i, where Xj+i = Xi + 1, equation (|16D with the ansatz (|36D 
and the relation (pT]) give us the generalization of relation (|3^), namely 

a:;::!}::?::^:: = - e s^;^\kp^,kp^,M:::!^^k... Q.,Qm = o,i(38) 

Q[,Q'2=0 

with 5* given by (^). Inserting the ansatz ( P^D in the boundary condition 

/(xi, Qi;...;Xn, Qn) = fix2, Q2;---; Xn, Qn] Xi + L, Qi) (39) 
we obtain the additional relation 

^Ql,---,Qn _ ikp-^L j^Q2,--,Qn,Ql (4:0) 
Pl,---,Pn P2,---,Pn,Pl V / 

which together with (^8]) should give us the eigenenergies. 

Successive applications of (|38| ) give us in general distinct relations be- 
tween the amplitudes. They are consistent because as we will see below 
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the 5* matrix coincides with the famous 6-vertex R matrix and satisfies 
the Yang-Baxter equation. Hence we may use these relations to obtain the 
eigenenergies of the Hamiltonian (p!3|). Applying the relation ( pSf ) n times 
we obtain from ( ^0]) a relation connecting the amplitudes with the same 
momenta, namely, 

j^Qi,...,Qn _ ^ikp-^L j^Q2,---,Qn,Qi _ / ^\n—l^ikp-^L 

E E s^^ikp,,kp;)^{kp,,kp^... (41) 

Q'„...,% Q'/,...,Q;( 

where we have introduced the harmless extra sum 

1= E ^QlQ'M{,Q.= E S^^{kp„kpJ. (42) 

Q'{,Q2=0 Q'{,Q2=0 

In order to fix the values of {kj} we should solve (PD, i.e., we should find 
the eigenvalues A{k) of the matrix 

T(A:)|« = E \(J[S^^\kp,,k:^ sgrg(A:p„, fc)| (43) 

Q" Q''=0 I \«=1 / J 



and the Bethe-ansatz equations which fix the set {ki} will be given from (^Tf ) 
by 



e-^.^ = (_l)-iA(A:^.) j = l,...,n. (44) 

We identify T(fc) as the transfer matrix of an inhomogeneous 6-vertex model, 
on a periodic lattice, with Boltzmann weights SQ^'Qf^kp^, k) (/ = !,... ,n). 
Consequently, in order to obtain the eigenenergies of the quantum chain (p!3[ ) 
we should diagonalize the above transfer matrix T{k). 

It is convenient to introduce the variables {uj} as in the NBAE of §2: 

e-- = (45) 
sin(M + 7r/3) ^ ^ 



In terms of these variables the S-matrix (^3f) have a difference form 

- U2) = Slj{ui - U2) 



,0,0/ N _ ci.V.. ^, ^ _ sin(Mi - M2 + vr/3) 



sin('U2 — Ml + 7r/3) ' 
11 



y3ei(«i-n2) 
2 sin(M2 — Ui + 7r/3) 

^ei(«2-«i) 
2 sin(M2 — Ml + 7r/3) 



(46) 



sin(M2 — Ui + 7r/3) 

where we should remind that we are considering q = exp(2i7r/3). 
Introducing for convenience the new amphtudes 

n 

^::3: = E ^q.o «p.)^?v; J: (47) 



we obtain instead of (|38D a similar relation with A instead of A and with the 
S'-matrix of the symmetric 6-vertex model, whose non-zero components are 
given by 

Sofii^) = Sljiu) = p sin(7r/3 - u) 

'S'o,'i('") = Sl'o{u) = p sinu (48) 

-^1,0 (m) = Sof{u) = p simi/3 

where p = 1/ sin(7r/3 + u). 

The matrix (^3D is the transfer matrix of an inhomogeneous 6-vertex 
model on the square lattice of width n = no + ni. In the next section we 
show that for the case where hq = k and rii = k + 1 (sector (k,k-|-l) of the 
6-vertex model on the lattice of width 2k + 1) this transfer matrix has the 
special eigenvalue A = 1 independently on the values of the parameters u 
and Uj {j = 1, . . . ,2n + 1). Consequently, since n = 2k + 1, (^^ reduces to 

e-*'^^ = l, j = l,...,n, (49) 

and the associated energies of the quantum chain ( pTSj) are free-fermion-like. 



4. The special eigenvalue 



Applying the Bethe ansatz method to the auxiliary transfer matrix T{k) = 
T{u) introduced in (^) one obtain the well known NBAE 10|, which were 
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considered in our paper . We have shown there that these NBAE are con- 
sistent with the FFC (^91) for the sectors (no, ni,n2) = {k, k + 1, L — 2k — 1). 
However it is more convenient here to follow an early paper of Baxter 
who considered the most general integrable inhomogeneous 6-vertex model. 
Any eigenvalue T{u) of the inhomogeneous model with p = 1 in (^) satisfy 
the equation 

T{u) Q{u) = (f[ sin(7r/3 -u + uj)] Q{u - 27r/3) 



+ 1 n sin(M - Uj) ) Q{u + 27r/3) (50) 
where Q{u) is an auxiliary trigonometric polynomial of degree tt-q, namely, 

no 

Q{u) = \\_sm{u — aj). (51) 



It is clear from (|50|) that T{u) is a trigonometric polynomial of degree n and 
from (^8]) the eigenvalues of (H3]) (where now p 7^ 1) are given by 



K{u) =T{u)/\{sin{'K/'i + u-Uj). (52) 
i=i 

This last expression implies that we have an eigenvalue A('u) = 1 if 

n 

T{u) = Y[ sin(7r/3 + u — Uj). (53) 

i=i 

The existence of this special eigenvalue was argued by Baxter for the more 
general case of the 8- vertex model with special values of the crossing param- 
eters [0. We can prove the existence of this eigenvalue for lattices with odd 
values of its width n. In this case we can rewrite Baxter's T-Q equation (50) 
as 

f{u) + f{u + 27r/3) + f{u - 2n/3) = (54) 

where 

n 

f{u) = Q{u + 27r/3) n sin(M - Uj). (55) 
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It is now clear that f{u) is a trigonometric polynomial of the degree n + nQ. 
Equation ( ^5]) coincides with (6) in where f{u) is also a trigonometric 
polynomial. 

In m it is shown that for any set of complex numbers Uj, j = 1, . . . ,2k + l, 
there exists a trigonometric polynomial Z[u) of degree n on the variable u 
such that 

2k+l 

f{u) = Z{u) sin(M — Uj) (56) 
i=i 

satisfies (|5^) p|. The degree of / is equal to n + no = 3/c + 1 and the degree of 
Z is equal to the degree of Q, so that k = uq. We obtain consequently that 
the construction of paper correspond to uq = k and n = no + ni = 2k + l, 
i.e. ni = k + 1, that complete the proof of existence of the special eigenvalue 



([53|) for odd values of n. 



5. The special wave function of the inhomogeneous 6-vertex 

model at g^JTr/s 

We consider several families of transfer matrices (|43|) corresponding to 
each distinct permutation P = {Pi,...,P„}. It follows from that 
Ap^''"'p" are the 2" components (Qi, ■■■,Qn = 0, 1) of an eigenvector of the 
transfer matrix (^31). Now we are going to investigate these components for 
our special eigenvalue A = 1 using the generalization arbitrary number of 
particles n, i. e., 

p={Pi,Pi+i} 

Y: crp,,p,+i ^: = g = i,2 (57) 

p={Pi,Pi+i} 

p={Pi,Pi+i} 

where all indices shown by dots are fixed. When we constructed the S'-matrix 
in ( P^D we expressed A'^^'p^ as a linear combination of Ap^'p^ and A^^'p^. 

'^Z{u) is the partition function of the inhomogeneous 6-vertex model with domain wall 
boundary conditions and with rapidities {u, Uj}, j = 1, . . . , 2fc + 1. 
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Now, on the other hand, we intend to express A'p^'p^ as a hnear combination 
of A^'^'^^^ and Ap^'p^ {Qi ^ Q2). Combining the first and the third equation 
of the set ( pTj) we obtain 

^p°P2 (e'^''^ - e'^'^i ) = (1 - qe'^P^ ) (1 - q-^d''^^ )^f,,p^ 

+ ^P.M%,P, (58) 

I /1 1.0 

+ ^P2,Pl^P2,Pl- 

Changing the variables as in (531) these equations are replaced by 



+ sin(Mp, - + 7r/3)A^2,Pi (^9) 

^Pl,P2^^^y^P2 ^Pl) — ^Pl,P2 

+ sin(Mpi -up,^+ 7r/3)y4^2°Pi- 



Generalizing these equations for arbitrary n and considering A instead of A 
(as in (^71) ) we obtain 

;i--,o,i,- 



v4"''p°p" sinfup,^, — ^^p,) = -^A' p p 

■••iPi,P;+iv ^ ^'+1 ^i'' 2 



+ sin(Mp, - Mp,^^ + 7r/3)Ajp^^^'p^^ = 



^:;p„p;;i,...sin(Mp,+, -Mpj = — a::p,>;^,,... (eo) 

+ sin(up, - Mp,_^^ + 7r/3)i:::;J^°;^;p,_,„ = 0. 
We have also a generalization of the second equation in (|27|) 



sinK^, -MP, +7r/3) A-^;'^;;;^ = 

- sin(«p, - up,^, + 7r/3) i^i^/J.p,,... (Q = 0, 1). (61) 
Moreover due to FFC (^) equation ( ^0]) leads to the cyclic symmetry 

AQl<---<Qn AQ2,---,Qn,Ql f RO^ 

^Pl,...,P„ - ^P2,...,Pn,Pl 
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which is also vahd for A. 

Let us begin with the simplest nontrivial case k = 1, where the sector 
of appearance of free-fermion levels is (1,2,L — 3). Equation (|T^ for the 
eigenvectors can be written as follows 

|no, rii >= |1, 2 >= 

J2 J2 fixi,Qi;x2,Q2;x3,Q3)\xi,Qi;x2,Q2;x3,Q3> (63) 

{Q} l<xi,X2,X3<L 

where we sum over the three sequences of {Q}: {0, 1, 1}, {1, 0, 1} and {1, 1, 0} 
that correspond to the configurations where Uq = 1, rii = 2 {n = rio+rii = 3). 
The amplitudes are given by the ansatz (|36|) 

fix,,Q,;X2,Q2;xs,Q3) = ^A«^;g^^^^e^(^^i^i+^^2-2+fcp3-3^ (64) 

p 

In the above expression there are 3 types of parameters {A} which are related 
among themselves by the cyclic symmetry: 

70,1,1 _ 71,1,0 _ 71,0,1 (r^.^ 

Before proceding let us introduce the simplified notations 

sin(Mp^ — up^ + 7r/3) 



■5Pl,P2 



sin(7r/3) 



sm(np, - -»pj 
sm(7r/3) 



The relation ( [61| ) can be reduced to the equation 

^0.1,1 _ 40,1,1 /a^x 

SP'i,P2 ^Pl,P2,P3 ~ ~^P2,P-i ^Pl,P3,P2 ) 



SO that due to (|55|) 

iS],,,P3 = ±C{P^}sp,,p, (68) 

where the sign depends on the parity of the permutation P = {Pi, P2, P3} 
and C{i},i = 1,2,3 are unknown coefficients. 

The relations (^) give us in particular the equation 

^Pi°,P2,P3'^'P2,Pl = ^Pl^,P2,P3 + ^^'l:^'2^P2^,Pl,P3- (^^) 
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Using the cyclic symmetry and ( |55D we obtain 

C{P2}sp,,p,dp2^P, = C{Pi}sp2,P3 - C{P2}sPi,P2Sp,,P3 (70) 
that together with the identity 

sp3,Pidp2,Pi + sp^,P2Sp^,P3 = sp^^p^ (71) 
give us C{P2} = C{P2} = C, i.e., up to a normahzation factor we have 

iSi,P3 = ±^P2,P3 (72) 

where the sign depends on the parity of the permutation P. 

Consider further the next sector (2,3,L — 5), i.e., uq = 2, rii = 3 and 
n = no + rii = 5. In this case we have two sets of {^4}, which are related due 
to the cyclic symmetry, namely. 



20,0,1,1,1 _ 70,1,1,1,0 _ _ 71,0,0,1,1 

"^Pl,P2,P3,P4,P5 ~ ^P2,P3,Pa,Ps,Pi " • • • " ^Ps ,Pl ,^2 ,^3,^4 

70,1,0,1,1 _ 71,0,1,1,0 _ _ 71,0,1,0,1 /7Q\ 

^Pl,P2,P3,P4,P5 ~ ^P2,P3,P4,P5,Pl " • • • " ^Ps ,Pl ,P2 ,P3,P4 " \''^) 

Let us begin with the first set. The relation (^) is solved by the ansatz 

^P°,P;!p3,P4,P5 = ±<^{^l'^2}SPi,P2Sp3,p,Sp3,P5Sp,,Pg (74) 

where the sign depends on the parity of the permutation P = {Pi, P2, P3, P^, P5} 
and C{i, j},i, j = 1,2,3,4,5 are symmetric unknown coefficients C{i,j} = 
C{j, i}. The system ( pDD contains in particular the three following equations: 

70,1,0,1,1 , _ 70,0,1,1,1 , ^ /iO,0,l,l,l 

"^Pl,P2,P3,P4,P5"'P3,P2 — ^Pl,P2,P3,P4,P5 ■^P2,P3^Pl,P3,P2,P4,P5 

70,1,0.1,1 7 _ 70,0,1,1,1 I /lO,0,l,l,l f^r\ 

^P2,Pl,P3,P4,P5^'P3,Pl — ^P2,Pl,P3,P4,P5 + '^Pl .Ps ^P2 ,P3,Pl ,P4,P5 V'^'' 

71,0,0,1,1 7 _ 70,1,0,1,1 I 70,1,0,1,1 

^Pl,P2,P3,P4,P5^'P2,Pl — ^Pl,P2,P3,P4,P5 + '^Pl ,P2 ^P2 ,Pl ,P3 ,P4,P5 " 

Excluding from this system ^4°'^'°'^'^, using the cyclic symmetry 

71,0,0,1,1 _ 70,0,1,1,1 

'^Pl,P2,P3,P4,P5 ~ '^P2,P3,P4,P5,Pl 

and limiting ourselves with the unit permutation we get 

70,0,1,1.1 J J J _ ( 70,0,1,1.1 I 70.0,1, 1,1\ 7 I 

^2,3,4,5,l"2,l«3,2«3,l — 1^1,2,3,4,5 + ■S2,3^1,3,2,4,5J"3,1 + 

/ 70,0,1,1,1 I 70,0.1,1,1\ 7 /^n\ 

1^2,1,3,4,5 + -51,3^2,3, l,4,5j"3,2'Sl, 2- V^) 
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Inserting here the ansatz (^) and using the antisymmetry of dij we obtain 

C{l,2}si,2S3,4S3,5(c?2,3S2,l - C^l.s) + C{l,3}di^3Si^3S2,3S2,4.S2,5 + 
C{2, 3}ci2,3S2,3('^l,2C?l,3-S4,lS5,l " Si,2Si,3Si,4Si,5) = 0. (77) 

Using the identity 

C^2,3S2,l - di^3 = -rfi,2S2,3 (78) 



and removing a common multipher ^2,3 from (|77D we obtain the more simple 
equation 

—C{1, 2}(ii 2Si^2S3,4S3,5 + C{1, 3}(ii_3Si 3S2,4S2,5 + 

C{2, 3}(i2,3(rfl,2rfl,3S4,l'S5,l - Sl,2Si,3Si,4Si,5) = 0. (79) 

By interchaging indices 1 and 2 we obtain a distinct equation. If we now 
exclude C{1, 2} from these two equations we get a relation between C{2, 3} 
and C{1,3}, namely, 

C{2, 3}c?2,3((ii^2C^l,3'S2,lS4,lS5,i + Si^2'5l,4'5l,5(s2,3 ~ ■52,1-51,3)) = 

C{1, 3}(ii,3((ii,2(i2,3Sl,2S4,2S5,2 + ^2,1 52,4^2,5 (^1,252,3 " Si^s)). (80) 

This last relation takes a nice form if we use the identities 

S2,3 - 52,151,3 = ^1,2(^1,3 -31,252,3 - Si,3 = 6^1,2^2,3 (81) 

and remove the common factors di 2, c?i.3 and (^2,3, i. e., 

C{2,3}(S2,1-S4,1-S5,1 + -Sl,2-Sl,4-Sl,5) = 

C{1, 3}(Si,2.S4,2-S5,2 + -52,1-32,4-52,5)- (82) 

Using standard trigonometric identities one can show that the left side com- 
bination 

S2,lS4,lS5,l + Si,2Si,4Si,5 = -|cOS(Ui + U2 - - U5) + 

cos(mi — U2 + U4 — M5) + cos(ui — U2 — U4 + U5)} (83) 
has a 5*4 symmetry, and consequently (^2|) reduces to 

C{2,3} = C{1,3}. (84) 
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This means that C does not depend on its indices and we have up to a 
normahzation factor 

^Pi°P2!p3,P4,P5 = ='='Spi,P2"5P3,P4"5P3,P5"5P4,P5- {^^) 

The second set of amphtudes can be found from the first relation in (ffSl). 
For example 

70,1,0,1.1 _ 7-1 / N _ 

^1,2,3,4,5 ~ "2,3"54,5l'5l,3'52,3"52,4'S2,5 ~ 'Sl,2'S3,4"S3,5 j = 
^^2i3'54,5{s2,3S2,5(si,3S2,4 - Si,2S3,4) + Si,2S3,4(s2,3S2,5 - ^3,5)} (86) 

and using the identities 

Sl,3S2,4 - Sl,2S3,4 = (^2,354,1 52,3^2,5 " ^3,5 = C?2,3S5,2 (87) 

we obtain 

^l!2;3A5 = 'S4,5(S2,3S2,5S4,1 + Si,2S3,4S5,2) , (88) 

or equivalently by using some additional identities we get 

^i^2fiAl = -54,5(52,352,455,1 + 51,253,554,2). (89) 
In a similar way we can derive the general answer 

70,1,0.1,1 / _ I _ 

^Pl,P2,P3,P4,P5/ ^Pi'Ps — 5P2,P3 5P2,P5 5p4,Pi + S p^^p^S p^^p^S ^p^ — 

5P2,P3 5P2,P4 5P5,Pi + 5Pi,P25p3,PgSp4,P2. (90) 

Equations (0) and ( pS] ) induce us to conjecture that for an arbitrary 
n the amplitudes {A} of the special wave function of the inhomogeneous 
6- vertex model are given by the ansatz 

^p;:.'.°,p:a+i,...,p2.+i = n ^p^,pJ n ^p.,p, (91) 

^<i<j<k k+l<i<j<2k+l 

where k and k + 1 are the numbers of particles of species and 1, respectively. 
We checked this formula for n = 7 analytically and for n = 9 using a brute- 
force diagonalization. It is a challenge to prove the validity of this formula 
for an arbitrary odd number n. The remaining amplitudes can be found by 
using this ansatz and sucessive application of the "recursion" relation (|60D. 
This completes our discussion of the SU{2) periodic case. 
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6. Summary and Conclusions 



In previous papers |1|, it was shown the existence of free-fermion-like 
energies for the anisotropic SU{N) Perk-Schultz model with anisotropy pa- 
rameter q = exp i7c{N — 1)/N. These solutions were found for general values 
of in the case of free boundary condition, where the model is SUq{N) 
invariant and for the SU{3) case in the periodic case. 

In §2 of this paper we show that the above observations, for the case 
of free boundaries, could be explained by a conjecture stating the existence 
of a special factorizable eigenvalue of the auxiliary inhomogeneous transfer 
matrix of a SUq{N — 1) vertex model with the same value of the anisotropy. 
Although we believe that a general derivation of such factorizable eigenvalue 
would be possible we restricted our analytical work in the simplest case of 
the periodic SU{3) Perk-Schultz model at g = expi^. In this case the 
associated transfer matrix is that of the inhomogeneous 6-vertex model and 
the existence of the factorizable eigenvalue, as shown in §4, follows from the 
T-Q Baxter equation (|50D. 

In §3 we review the coordinate Bethe ansatz and show how to relate the 
wave function components of the eigenvectors of the SU{3) quantum chain 
with periodic boundaries in terms of the components of the eigenvectors of 
the inhomogeneous 6-vertex model. In particular all the free-fermion-like 
energies are related to a single factorizable eigenvalue of the inhomogeneous 
6-vertex model. 

In §5 exploring the existence of the free-fermion-like solutions of the SU{3) 
chain at g = expi^ we show how to produce the recurrence relations that 
allows the computation of the wave vector amplitudes related to the special 
factorizable eigenvector of the 6-vertex model. These relations, although not 
simple, give us a systematic way to derive all the eigenf unctions of the free- 
fermion part of the eigenspectrum of the periodic SU (3) Perk-Schultz model 
at g = exp i^. 

As an application let us consider the free-fermion branch of this last 
model in the sector Sk = {k, k + 1, L — 2k — 1). From (JT^) the corresponding 
eigenenergies are given by 

E, = -^(l + 2cos(27rj7L)), 

where / is any subset of (1,2,..., L) with n = 2k + 1 distinct elements. 
Enumerating these elements by the index a = 1,2, . . . ,n, the wave function 
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is given by 



fix,, Qi; . . . Q„) = Y: J;e^('=-i-^+---+'=-"-"), 

p 



where ka = 2iTja/L, a = (1, 2, . . . , n) are the momenta of the elementary 
free-fermion excitations. 

Let us hmit ourselves to the subsets / with elements ja, satisfying the 
constraint j < m oi L — m < j, where m is a positive integer. Due to 
conjecture 2 of (0) the lowest eigenenergy in the sector Sk = {k, k + 1, L — 
2k — 1) belongs to this part of spectrum if we choose m > k. 

Now we fix k and m with m > k and consider the bulk limit L — > oo. Due 
to above mentioned constrains for the values of k we have two possibilities: 
ka ^ or ka ^ 271. Consequently from (^) all the parameters Ua that 
fix the auxiliary 6- vertex model become equal to 7r/2 and we obtain the 
homogeneous model with the special eigenvector A^^''"'^" which (up to a 
sign factor) does not depend on the particular permutation P\ The wave 
function can then be written as 

Qi; • • • ; Vn, Qn) = A^^-"^- ^(_i)Pe2-(^-i^i+-+^--"^"), (92) 

p 

where we introduced the new coordinates yi = Xi/L {i = 1, . . . ,n). In partic- 
ular this result shows that in the sectors {k, k + 1, oo) there exist a lot of eigen- 
states (including the one with lowest eigenenergy in the sector) whose wave 
functions components are given by the product of Slater determinants and the 
components of the ground state wave function of the XXZ model with 2A; + 1 
sites. Let us consider for example the sector ^2 = (2, 3, L — 5). In order to 
obtain the lowest eigenenergy in this sector we chose / = 1, 2; L — 2, L — 1, L. 
The Slater determinant in ( p^ ) reduces to the Vandermonde determinant and 
we obtain (up to a normalization) 

fiy,,Q,;...,y,,Q,)= H sin 7r(|/, - |/fc)A«^■■••''3^ (93) 

l<j<fe<5 

where we have 

^00111 _ ^01110 _ ^11100 _ ^11001 _ ^10011 _ 
^01011 _ ^10110 _ ^01101 _ ^11010 _ ^10101 _ 2 

From these components we see, for example, that the probability to find 
0-particles , separated by 1-particles is equal to 2 * 2/(1 * 1 + 2 * 2) =4/5. 
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Moreover, it is important mention that when the inhomogeneity of the 
auxihary 6- vertex model disappears, and the wave function that corresponds 
to the special factorizable eigenvalue is the same as that of the ground state 
of the XXZ spin chain with A = —1/2, which possess for L = 2r;, + 1 quite 
interesting combinatorial properties. 

Due to the conjecture announced at the end of section 2, we have a 
generalization of this special wave function to the SU{N — 1) - invariant 
case {q = e*^^"^)'^/^). So we suspect that the ground state function of the 
SU{N — 1) - invariant spin chain can also exhibit interesting combinatorial 
properties for the special value q = e*(^~^)'^/^. Indeed, using a bruteforce 
diagonalization of these quantum invariant chains we have found that the 
ratio defined by 

R, = (94) 

where {vi} are the wave function components of the ground state, has a 
simple form depending on the boundary condition and on the parity of the 
lattice size. For the L-sites XXZ spin chain at A = —1/2 we have Rl = \/3 
0, where a = L — 1 or a = L depending if the length L of the chain is odd 
or even, respectively. In the case where L is odd the chain has the boundary 
condition periodic or SUq{2) invariant, and for even values of L the chain 
has boundary condition of twisted type or a SUq{2) invariant one. We can 
present these results in a compact form: 



Rl+2IRl = 3, Ri = 1, R2 = 3. (95) 

The numerical results coming from bruteforce numerical diagonalizations 
of chains with SUq{N) symmetry (free boundary condition) followed by a 
fitting with special irrational numbers give us the values of Rl for the next 
cases: 

• SUg{3) with q = -e^^l^ 2<L<9 

Rl+3/Rl = (1 + v^)' (96) 
i?i = 1 R2 = l + V2 i?3 = (1 + V2f- 
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• SUq{A) with q = -e'^l^ 2<L<8 

Rl+a/Rl = (5 + 2^5)2 (97) 
R^ = l R2 = V5 i?3 = 5 + 2^5 i?4 = (5 + 2v^)l 

• SUg{5) and q = -e^^l^ 2<L<7 

Rl+^/Rl = (2 + v^)^ 

Ri = l i?2 = (2 + v^)/V3 i?3 = (2 + V3)VV3 (98) 
/24 = (2 + VS)^ i?5 = (2 + VS)^ 

The numbers Rl, L < N which are necessary for the use of the recursion 
relations can be found from the exphcit expression of the corresponding wave 
function. Using an approach, described in Appendix B of paper one can 
show that for L < N: 

«.-(l^)"'nl^ (99) 

For the cases where L > N although we cannot prove, our numerical results 
indicate the nice recursion relation 

Rl+n-1 = Rn-1 Rl- (100) 
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